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ABSTRACT 

We investigate the dependence of the galaxy infrared luminosity function (LF) and the associated 
stellar mass function (SMF) on environment and spectral type using photometry from the Two Micron 
All Sky Survey and redshifts from the Las Campanas Redshift Survey for galaxies brighter than Mj < 
—19 + 51og/i. In the field environment, galaxies with emission lines have LFs with much steeper faint 
end slopes (aj = —1.39) than galaxies without emission lines (a,/ = —0.59). In the cluster environment, 
however, even the non-emission line galaxies have a steep faint-end LF (aj = —1.22). There is also a 
significant (95%) difference between the overall cluster and field LFs, Aaj = -0.34, AM] = —0.54. All 
of these variations are more pronounced in the SMFs, which we compute by relating the strength of the 
4000A break in the optical spectra to a stellar mass-to-light ratio. 

Subject headings: galaxies:clusters:general — galaxiesduminosity function — galaxies:mass 
function — galaxies : evolution — infrared : galaxies 



1. INTRODUCTION 



Boselli [1995; Bel1 & de Jon S US EH et al 



The shape of the stellar mass function (SMF) ™ID and f ne f rl y ^sensitive to the presence of 



is a fundamental property of the galaxy popula- 
tion because it constrains the baryon content of 
galaxies, the fraction of all baryons in stars, and 
the chemical evolution of galaxies. Until recently, 
however, observational work has been unable to 
provide a reliable, direct estimate of the SMF. 
This failure arises primarily because the conver- 
sion of the luminosity function (LF) into the SMF 
requires the poorly-constrained optical mass-to- 
light (M/L) ratio, which in turn depends sen- 
sitively on both the stellar populations and the 
spatial distribution of dust and stars, neither of 
which is known for a substantial sample of galax- 
ies. Complex models of dust and star formation 
are therefore required to compare the observations 
with theoretical pre dictions ( |Cole et al. 2000a| ; 



Granato et al, 2000| ), 



dust at the levels expected in most galaxies (the 
typical internal extinction at K s is ~ 0.05 mag; 
|Silva et al. 19981 ). 

Three recent developments motivate us to at- 
tempt the measurement of the SMF's environ- 
mental dependence. First, large IR surveys are 
an entirely new resource, due to the advent of 
large format detectors, dedicated survey tele- 
scopes, and significant human and material re- 
sources. Second, the availability and uniform 
analysis of large optical and IR surveys has led 
to a consensus on the shape of the global galaxy 
LF at both optical and IR wavelengths. Although 
controversy persisted for many years regarding 
differences between the optical LFs obtained by 
different investigators ( e.g., |Loveday et al. 1992j; 
Efstathiou et al. 1988); |Marzke et al. 1994b|; |Lin 



et al. 1996| ;~|Z ucca et al. 1997] ; [Ratcliffe et al 
1998|) , these differences have now largely been re- 



It is now possible to mitigate these difficulties 



by combining spectroscopic catalogs with large 



infrared (IR) surveys. The sensitivity of M/L 
to stellar populations and extinction, which has 
plagued the transformation of the optical LF to a 



solved (Blanton et a/. [2000| ). At IR wavelengths 
agreement among different samples has been ex- 
traordinarily good (IMobasher et al. 1993; IGlaze- 



SMF, is lessened at IR wavelengths. M/L varies 



brook et al. 1995|; |Gardner et al. 1997 ; gzokoly 



et al. 1998|; |Loveday 2000| ; |Kochanek et al. 2000b 
Cole et al. 2000b| ). This long-awaited concordance 



by only a factor of ^ two over a large range of 



star formation histories 



v e.g., Gavazzi, Pierini & 



1 This publication makes use of data products form the Two Micron All Sky Survey (2MASS), which is a joint project of the University of 
Massachusetts and the Infrared Processing and Analysis Center /California Institute of Technology, funded by the National Aeronautics and 
Space Administration and the National Science Foundation. 
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makes it possible to confidently proceed and in- 
vestigate variations in the LF for different galaxy 
types, different environments, and different red- 
shifts. Third, the availability of large, cosmologi- 
cal simulations allow the dependence of the SMF 
on environment to be modeled with some reliabil- 
ity. However, the translation from the dark mat- 
ter halo mass function to the SMF still depends on 
such complicated and ill-understood factors as the 
bias function ( |White et al. 1987); [Narayanan et aT. 
20001} , the gas cooling rate (|Rees and Ostriker 
1977|), and reheating due to feedback from star 
formation (e.g., Larson |1974| ; Cole |1991| ; Balogh 
et aZ. pOOlp . 

A dependence of the SMF on galaxy type would 
not be surprising, due to the distinctly different 
star formation histories of late- and early-type 
galaxies. In addition, we can expect an environ- 
mental dependence of the SMF because (1) the 
merging history (e.g., Lacey and Cole 1993|) and 
bias level ( |Bardeen et al. 1986| ; [Narayanan et aT 



20000 of halos in clusters are different from those 
of isolated halos, (2) there are physical processes, 
such as galaxy harassment (Moore et al. 1999|) , 
that may operate in dense cluster environments 
and alter the halo mass function, and (3) the 
stellar populations of cluster galaxies are differ- 
ent from those of field galaxies, perhaps indicative 
of different star formation histories (e.g., |Dresslei 
et all 19851; [Balogh et al. 19971; Balogh et al 



Poggianti et al. 1999j ; [Moss and Whittle 



1998 
2000 

Such expectations have already received some 
support from observations. It is well known that 
the optical LF depends on galaxy morphology 
e.g., [Loveday et al 1992|; Binggcli et al 1988| ; 



Marzke et al. 1994a| [Marzke et al. 1998|) , pri- 
marily in the sense that it steepens toward later 
galaxy types. Analogous trends are also seen if 
the population is divided by emission line strength 
([Ellis et al. 1996|; |Lin et al. 19"96|; |Lovcday 



et al 



1998 



1999|; IChristlein 200q), color (Marzke anc 



da C osta 1997| [Lin et al. 1996| ; |Trcyer et al 



or principal spectral components ( Brom 



ley e > al. 1998j ; |Folkes et al. Tggg). In addition 



to the dependence on type, optical observations 
also indicate that the dwarf-to-giant ratio steadily 
increases with increasing environmental density 
([Sandage et al. 1985; |de Propris et al. 1995|; [Fer-j 



ffusoih and Sandage 1991|; Bromley et al. 19981; 



Zabludoft and Mulchaey 2000| ; IChristlein 2000| ). 

The trends seen in the optical LFs must be ex- 
plored in the IR to determine whether they reflect 
underlying variations in the SMF or in M/L. In 



one recent example of such work, Kochanek et 
al. ( [2000 b| ) find that the K s — band LF of galaxies 
has the same dependence on morphology as seen 
at optical wavelengths. This finding is strong evi- 
dence that the galactic SMF is not universal, and 
that early-type galaxies cannot simply be an un- 
biased subsample of late-type galaxies in which 
star formation has ceased. Our primary purpose 
is to further explore in the IR the trends found in 
the optical. The compilation of the data we use 
is described in §§. In §^ we calculate the J and 
K s band galaxy luminosity functions ClS db func- 
tion of environment and spectral type. We find 
that the trends in optical data are also present in 
the IR data, and more significantly in the stellar 
mass function. The implications are explored in 
§f|, and we summarize our conclusions in §|5|. 

2. DATA 

To construct infrared luminosity functions, we 
use photometry from the Two Micron All Sky Sur- 
vey (2MASS; |Jarrett et al. 2000|) , and redshifts 
from the Las Campanas Redshift Survey (LCRS; 
Shectman et al. 1996|) . The LCRS consists of over 



25,000 galaxies with redshifts, approximately cov- 
ering the magnitude range R c = 15 to R c = 17.7, 
where R c is the "hybrid" Kron- Cousins R system 
described in Lin et al. ( |1996| ). The survey in- 
cludes some fields observed with a 50 fiber config- 
uration, and other, later fields observed with a 112 
fiber configuration. The 50 fiber data are prob- 
lematic because the sampling fraction is lower 
and the apparent magnitude range of the tar- 
geted galaxies is typically less than 1.5 magni- 
tudes. Galaxies with high central surface bright- 
nesses are excluded from the redshift survey, and 
this limit is about 0.7 mag fainter in the 50 fiber 
data. To ensure a homogeneous and maximally 
complete sample, we restrict ourselves to the 112 
fiber data, as did Lin et al. ( |1996| ) for their fiducial 
computation of the optical LF. 

We use the second incremental data release of 
the 2MASS catalog, which contains photometry 
in J, H and K s bands of both extended and point 
sources. The overlap of this release and the 112- 
fiber LCRS fields is patchy, covering only about 
40% of the LCRS. Cross-correlation between the 
catalogs is done by matching galaxy astromet- 
ric positions. Following the example of Cole et 
al. ( P000b| ), we allow matches between galaxies 
whose positions differ by less than 0.75 times the 
J-band Kron radius, to allow for the greater po- 
sitional uncertainty of larger galaxies. 
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2.1. Spectral Classification 

Zabludoff et al. ( |1996| ) measured the rest frame 
equivalent width of the [OII]A3727A emission line, 
W^OII), and the strength of the Balmer break, 
-Dztooo- We use H^OII) to divide the sample into 
broad spectral classes: emission line (EL) galaxies 
with (Oil) > 5 A and non-emission line (NEL) 
galaxies with W o (0II)< 5A. Two caveats accom- 
pany this division. First, the [Oil] line correlates 
with Ha emission, albeit with considerable scat- 
ter, suggesting that it is a good indicator of star 
formation in an average sense (e.g., [Kennicutt 
19920. However, H^OII) also depends on metal- 
licity and ionization, and on extinction if dust is 
non-uniformly distributed (i.e., more prevalent in 
HII regions; |Charlot and Longhetti 20QT|) . All of 
these properties are expected to depend on galaxy 
luminosity. The correlation between [Oil] and Ha 
does indeed have a luminosity dependence ( [lresse| 
et al [ 1999] ; |Jansen et al. 2001| ). Second, al- 
though, the fibers used in the LCRS are quite 



large (3'.'5), Kochanek (|2000a|) suggests that aper- 
ture effects may still be important. We address 
this issue in 83. 1J. 



2.2. Environment Classification 

Group and cluster catalogs have been gener- 
ated by Christlein (|2000|) from the LCRS using a 
friends-of-friends algorithm to search for overden- 
sities in redshift space. For each association, the 
one dimensional velocity dispersion, oi, is com- 
puted. We define galaxy clusters to be those with 
0i > 400 km/s, corresponding to a virial temper- 
ature of about 1 keV. Associations with smaller 
Ox are termed "groups" and any galaxy not found 
to be in an overdense region is assigned to the 
"field" sample. Only associations with at least 
three members are retained in the final catalog; 
we discuss the possibility of biases related to this 
effect both below and in § [2.4| . 

Velocity dispersions computed from small sam- 
ples are biased low due to insufficient sampling 
of the wings of the distribution and the poor, bi- 
ased determina tion of the mean ( ^abludotf and 
Mulchaey 1998|) . This bias results in some frac- 
tion of the associations that properly belong in 
the "cluster" sample being incorrectly placed into 
the "group" sample. Furthermore, identification 
of groups with low velocity dispersions is sensitive 
to the linking length used in the cluster-finding 
algorithm, and such groups are more likely to be 
contaminated by interlopers. This contamination 
results in some groups that are false detections, 



and the galaxies in those "groups" should be prop- 
erly identified as isolated galaxies. To the degree 
that the environments are misclassified, any un- 
derlying differences among environments will be 
diluted. 

2.3. Completeness 

We compute the luminosity functions in both 
the K s band, which is commonly used, and the 
J band, where the 2MASS photometry is deep- 
est. In J, we use the Kron aperture magnitudes 
as the best estimate of the total magnitude. Cole 
et al. ( P000b| ) demonstrate that the most robust 
K s magnitudes are obtained by combining the 
Kron J-magnitudes with the default J — K s col- 
ors. This nonintuitive result arises because the 
J images have a higher signal-to-noise ratio and 
allow a more accurate determination of the Kron 
aperture size. We adopt this magnitude definition 
of K s . A small fraction (2.2%) of the galaxies in 
the LCRS within our magnitude range (defined 
below) are misclassified as stars in 2MASS. For 
these galaxies, we use the 8" aperture magnitudes 
in the point source catalog. 

It is important to check the overall complete- 
ness of the matched LCRS-2MASS catalog used 
in our analysis. In Figure [I| we show the fraction 
of the galaxies in 2MASS for which a successful 
match is made with a galaxy in the LCRS. We 
restrict our comparison to the largest (20° x 1.5°) 
contiguous patch covered by both surveys, which 
includes about 1300 galaxies. The LCRS is not 
intended to be complete, but is a random sample 
of a magnitude-limited survey. The dashed line 
shows the effect of applying the weights necessary 
to correct for this incompleteness. At K s > 12.2 
and J > 13.2, the corrected completeness is close 
to 100%. At brighter magnitudes, the complete- 
ness drops, due primarily to the LCRS's upper 
magnitude cutoff (R ~ 15); we therefore limit 
the sample to magnitudes fainter than this limit. 
Cole et al. ( |2000b|) show that 2MASS is highly 
complete and is not missing a substantial number 
of low surface brightness galaxies. By association 
there can be no significant incompleteness in low 
surface brightness galaxies in the matched LCRS- 
2MASS catalog, at least to the relatively bright 
lower magnitude limit of 2MASS. 

The 2MASS catalog is complete to K s = 13.2, 
and J = 14.5, as shown by Cole et al. ( |2000b| ). Up 
to about 0.5 mag fainter than these limits, the in- 
completeness is primarily due to misclassification 
of faint galaxies as stars and small biases in the 
magnitude measurements. We mitigate some of 
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Fig. 1. — The completeness of the LCRS, relative to 2MASS, as a function of J and K a magnitudes. The solid line shows the fraction 
of galaxies in the 2MASS catalog that have a corresponding match in the LCRS spectroscopic catalog. The dotted line shows the result of 
weighting the LCRS counts to account for the sampling fraction. The incompleteness at the bright end is due to the bright magnitude cutoff 
of the LCRS. 



this incompleteness by using the LCRS to iden- 
tify galaxies that 2MASS has classified as stars. 
Furthermore, because our sample is considerably 
smaller than that of Cole et al. , we are less con- 
cerned with systematic biasses of order ^0.1 mag, 
and we can maximize our useful sample size by 
considering fainter magnitude limits. 

One other way to test the completeness of 
our sample is by considering the distribution of 
^yVmax- For a galaxy at redshift z, V(z) is the 
volume between it and the lowest redshift z min at 
which it would still have been selected in the sam- 
ple (due to the bright magnitude limit). V milx is 
the volume between £ min and z max , where 2 max is 
the redshift at which the galaxy would drop out of 
the sample due to the faint magnitude limit. For 
a uniform spatial distribution of objects the mean 
value of V/V max is 0.5. Although galaxies are 
clustered, we expect that over the large volume 
surveyed here the distribution can be thought of 
as uniform. In Figure ^| we show the distribu- 
tion of V/V m3X for the K s band sample, limited 
at K s = 13.7, and the J band sample, limited 
at J = 15.0. With these limits, the distribu- 
tion of V/V m3X is approximately uniform, though 
there is a systematic variation of about 15%. The 
mean value of V/V ma , x is about 3a larger than the 
value of 0.5, which it would be for a statistically 
complete sampling of a spatially uniform popula- 
tion. This appears to be partly due to a deficit 



of galaxies with V/V max < 0.3, perhaps a conse- 
quence of Malmquist bias becoming important at 
K s > 13.5. Indeed, moving the magnitude limit 
brighter by 0.2 mag does make the distribution 
of V/V max more uniform. None of the trends we 
discuss in this paper are affected by decreasing 
the magnitude limit, although the uncertainties 
are increased and, consequently, the significance 
of the results is diminished. To maximize the use- 
ful sample size, we therefore restrict the K s sam- 
ple to 12.2 < K s < 13.7, which consists of 2673 
galaxies (including 274 cluster galaxies and 954 
group galaxies). In J, we consider the sample of 
3408 galaxies with 13.2 < J < 15.0 (including 346 
cluster galaxies and 1202 group galaxies). 

2.4. Potential Biases and Systematic Uncertainties 

Clusters and groups (associations) are defined 
by identifying friends-of-friends in redshift space. 
Only associations with three or more members 
are retained in the final catalog. A concern is 
that near this cutoff, associations with steeper 
faint end slopes a, or brighter characteristic mag- 
nitudes, might be preferentially included in the 
catalog because more member galaxies would be 
above the magnitude limit. This could bias the 
final LF; however, we expect the overall effect to 
be small because the majority of the associations 
in our catalogue have more than the minimum 
three members, and are unlikely to drop below 
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Fig. 2. — The distribution of V/V m a,x, for the K s band sample (solid line), limited at K s = 13.7, and the J band sample (dashed line) 
limited at J = 15.0. The mean value is 0.521±0.006 for the K s sample and 0.513 ±0.005 for the J sample. 



this threshold given reasonable variations in the 
LF. To test for this potential effect, we have con- 
structed mock catalogs, as described in Christlein 
Q2000D , by adding ~350,000 artificial associations 
to the LCRS catalog. In one catalog, the as- 
sociations have LFs with steep faint end slopes 
of a = —1.2, and in another they have shallow 
slopes, a = —0.6. As expected, we find no evi- 
dence from this experiment that the recovery frac- 
tion of associations with steeper faint-end slopes 
is higher than for those with flatter slopes. Thus, 
this potential bias does not have a significant ef- 
fect on our selection, even when considering mock 
samples containing many more galaxies than the 
actual data. 

There is a potential bias in any redshift cata- 
log to be more complete in emission line galaxies 
because it is easier to measure redshifts for these 
galaxies than for non-emission line galaxies. This 
bias is particularly acute for the faintest galax- 
ies in the sample. The magnitude ranges of the 
LCRS and 2MASS are not well matched, so only 
the brightest galaxies in the LCRS are included 
in our sample. In particular, most of the galaxies 
have 15 < R c < 16.5, well above the completeness 
limit of the LCRS. Because the magnitude range 
is small, and because we are far from the mag- 
nitude limit of the LCRS, there is little danger 
that a weak magnitude-dependent completeness 
function of the NEL galaxies has much effect on 
our analysis. In Figure || we show the V/V majX 



distributions for the NEL and EL galaxies in the 
K s selected sample. For the emission line galax- 
ies, (V/Vmax) = 0.515 ± 0.008, while for the NEL 
galaxies (V/V max ) = 0.525 ± 0.008. There is no 
evidence that the completeness of the EL sample 
differs from that of the NEL sample. 

3. RESULTS 

We compute the luminosity functions using 
both the non-parametric method of Efstathiou, 
Ellis & Peter son (|T988| , EEP) and by assuming 
a Schechter (|1976j ) function shape as described 
by Sandage, Tamman & Yahil ( |1979| , STY). Lu- 
minosities are calculated assuming an Q = 1, 
A = cosmology, with the Hubble constant 
parametrized as H = lOO/i km s _1 Mpc -1 . 
Where not explicitly shown, we use h = 1. K- 
corrections are taken from Poggianti ( |1997Q ; they 
are small (less than 0.1 mag at z = 0.1) and nearly 
independent of galaxy type. We neglect evolu- 
tionary corrections and extinction effects, both of 
which affect the magnitudes at less than the 0.1 
mag level (e.g., Cole et al. 2000b). 

We do not compute volume densities to nor- 
malize the luminosity functions. An estimate of 
the global normalization has already been pre- 
sented in Cole et al. ( |2000b| ). For the field galaxies 
and the total sample, we normalize the luminos- 
ity functions to the weighted number of galaxies 
brighter than M Ka = -21.5 and Mj = -20.5. For 
the group and cluster environments, we divide by 
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Fig. 3. — For the K a selected sample, the distribution of V/Vmax for emission line galaxies (solid line) and NEL galaxies (dashed line). 
The mean value is 0.515±0.008 for the emission line galaxies and 0.525±0.008 for the NEL galaxies. 



the number of associations to obtain the number 
of galaxies per group or cluster that are brighter 
than this limit. 

In Figure f| we show the J and i^-band lu- 
minosity functions for the full sample. Both 
are well fit by Schechter functions with param- 
eters M* K = -23.48 + 51og/i, a K = -1.10, and 
M*j= -22.23 + 51og/i, aj= -0.96. In the non- 
parametric method, the uncertainties are deter- 
mined from the information matrix (EEP). The 
uncertainties in the Schechter parameters are de- 
termined from contours of the likelihood function, 
and are shown as 67% confidence elliptical con- 
tours in Figure [5[ In Figure |^, we compare our 
result with recent measurements by Kochanek et 
al. ( pOOObl) and Cole et al. ( P000b|) . As shown in 
Cole et al. , the small difference between their re- 
sult and that of Kochanek et al. is mostly due 
to the difference in the magnitude definition (we 
use that of Cole et al.). For the K s sample, we 
are in good agreement with Cole et al. For the 
J sample, the la contours do not overlap, which 
may reflect the incompleteness introduced by in- 
cluding galaxies 0.5 magnitudes fainter than the 
nominal 2MASS completeness limit. We recover 
the Cole et al. result if we impose the same mag- 
nitude limits. While statistically significant, the 
difference in the J band results is small, and does 
not affect our conclusions based on subdividing 
the sample according to spectral type and envi- 
ronment. 



3.1. Dependence on Spectral Type 

In Figure ^] we present the LFs divided by spec- 
tral type. In J and K s the NEL galaxies have 
a significantly flatter faint-end slope than the EL 
galaxies; this result is also evident from the error 
ellipses in Fi gure |5| 

Kochanek ( |2000a|) has recently suggested that 
the fixed angular size of the fibers used to ob- 
tain the spectra leads to a substantial aperture 
related bias. To address the significance of this 
effect we show the fraction of NEL galaxies SIS db 
function of redshift, /nel(-z), in four bins of abso- 
lute magnitude (Figure 0). There is a weak trend, 
in the expected sense that more distant galaxies 
are more likely to have spectra with emission lines 
because the fixed aperture includes more of the 
galaxy (we have already shown, in § |2.4| , that this 
is not due to a redshift dependent incompleteness 
in NEL galaxies) . Although low luminosity galax- 
ies tend to be physically smaller than high lumi- 
nosity galaxies, they also lie at preferentially lower 
redshifts in a magnitude-limited sample. Thus, 
we expect the two effects to partially cancel out 
and reduce the size of this bias. 

3.2. Environmental Dependence 

As described in §Q, we divide our sample into 
three environmental categories, corresponding to 
field, groups, and clusters. In Figure [| we show 
the luminosity functions of these three classes of 
galaxies, compared with the best fit Schechter 
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Fig. 4. — The J (open symbols) and K s -band (filled symbols) luminosity functions. The points and solid lines are the determinations from 
the EEP and STY methods, respectively. The luminosity functions are arbitrarily normalized to the weighted number of galaxies brighter 
than Mk s = —21.5 and Mj = —20.5 (the J- band LF is arbitrarily shifted vertically for clarity). 
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Fig. 5. — The error ellipses (67% confidence contours) on the Schechter parameters of the luminosity functions, for the full sample (filled 
circle), EL galaxies (square) and NEL galaxies (cross). Recent determinations by Kochanek ct al. (2000) and Cole et al. (2000b) are marked 
K00 and COO respectively. The small difference between the COO and K00 results are due to the different magnitude systems used. All 
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Fig. 6. — The J and A" s -band luminosity functions divided by spectral type. Open squares represent emission line galaxies, while crosses 
represent non-emission line galaxies. The luminosity functions are arbitrarily normalized to the weighted number of galaxies brighter than 
M Ks = -21.5 and Mj = -20.5. 
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Fig. 7. — The fraction of NEL galaxies as a function of redshift is shown in the for four luminosity bins: —25.5 < Mjf s < —24.5 (filled 
circles); -24.5 < M Ks < -23.5 (open circles); -23.5 < M Kg < -22.5 (crosses); and -22.5 < M Ks < -21.5 (triangles). 
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functions for the full sample; the 67% confidence 
error ellipses on the Schechter parameters are 
shown in Figures and [10]. The Schechter pa- 



rameters of the J— band cluster LF are inconsis- 
tent with those of the field LF at the 95% con- 
fidence level, as determined from the likelihood 
contours. The trend is more significant in the 
J— band due to the larger sample size as a conse- 
quence of the deeper photometry. The difference 
is in the sense that the cluster LF has a brighter 
Mj and a steeper aj\ however, given the degener- 
acy between these parameters it is not clear pre- 
cisely how the LF shape varies. If we restrict 
the Schechter fit to a limited luminosity range 
(e.g., brighter than Mj = —20.5, or fainter than 
Mj = —23.5), the best- fit parameters change 
along the major axis of the error ellipse. How- 
ever, both the cluster and field parameters move 
in the same direction, and the difference between 
them is approximately maintained. It appears, 
therefore, that the environmental dependence of 
the LF is not restricted to just the bright or faint 
end. 

The differences among the LFs of different envi- 
ronments are larger when only the NEL galaxies 
are considered, as shown in Figure []. In the field, 
the LF of NEL galaxies is much flatter than that 
of the full sample, but in clusters the LFs of EL 
and NEL are statistically indistinguishable. We 
discuss this result further in §fO| The Schechter 
parameters of the J and K s luminosity functions, 
for each data subsample, are tabulated in Table 
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4. DISCUSSION 
4.1. Comparison with Literature 



It is encouraging that the global luminosity 
functions measured here, based on redshifts from 
the LCRS, are in such good agreement with re- 
cent determinations by Cole et al. ( |2000b|) and 
Kochanek et al. ( [2000b ). As shown by those au- 
thors, all infrared LFs derived from large samples 
are consistent (e.g., [CTobashcr et al. 1993|; plazc- 
brook et al. 1995|; Gardner et al. 1997| ; [Szokolyi 
et al. 1998; |Lovcday 20011 ). This concordance is 
in contrast with the historical situation at opti- 
cal wavelengths, where there has been consider- 
able disagreement between results in the litera- 
ture. In particular, optical LFs derived from the 
LCRS ( [Lin et al. 199B] |Christlein 20001 ) are gen- 
erally shallower than most other measurements, 
and it has recently been shown that this differ- 
ence results from an underestimation of galaxy 
magnitudes due to the shallow isophotal limit of 



the LCRS photometry and from the exclusion of 
low central surface brightness galaxies (Blanton et 
al. |2000[ ). The fact that we do not find a shallow 
slope in the infrared LF supports the interpre- 
tation that, at least over the relatively narrow, 
bright magnitude range of 2MASS, it is the R 
magnitudes, and not the survey selection itself, 
that is responsible for the discrepancy in the op- 
tical. 

Analysis of the optical LF dependence on en- 
vironment was done by Christlein ( [2000| ), using 
the LCRS sample. Although the LCRS magni- 
tudes are systematically biased due to the bright 
isophotal limit, this bias applies uniformly to the 
whole sample and is therefore unlikely to affect 
the trends discussed in that paper. Christlein 
found strong evolution of the LF with environ- 
ment, with the Schechter parameters changing 
along the major axis of the ellipse, toward steeper 
a and brighter M* with larger a. The magni- 
tude and sense of this change are fully consis- 
tent with the results presented in § |3.2| . The fact 
that Christlein finds such a continuous depen- 
dence on a (something we are unable to investi- 
gate due to our smaller sample size) provides en- 
couraging support that the effect is real. Zablud- 
off & Mulchaey (|2000| ) found similar evidence for 
steep, bright i?-band LFs in X-ray selected groups 
(a = —1.3). Because these groups generally have 
kT ^ 1 keV, many would be considered "clusters" 
by our criteria, and their results for steeper opti- 
cal LFs in these systems are very similar to our 
results for the cluster infrared LFs. 

Previously published evidence for an environ- 
mental dependence of the infrared LF is not as 
strong as for the optical LF. Both De Propris 
et al. (|T998D and Andreon & Pello (|2000D find 
a steep faint end slope in the H band LF of 
the Coma cluster, but only at magnitudes fainter 
than our sample probes. Trentham & Mobasher 
(|1998|) find no evidence for a difference between 
the i^-band LFs of cluster and field galaxies, over 
—24 < Mk < —22, although the uncertainties are 
large. Our sample is the first that is large enough 
to detect this weak effect at these luminosities. 
We look forward to dramatically improving the 
precision of this result when the full 2MASS cat- 
alog is released. 

4.2. The Stellar Mass Function 

Do the observed dependencies of the infrared 
LF on environment and spectral type necessarily 
imply an analogous difference in the stellar mass 
functions? The infrared light is a good, but not 
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Fig. 8. — The J and ii" s -band luminosity functions as a function of environment, shown as the error bars (EEP method) and the solid 
lines (STY method). The luminosity functions are arbitrarily normalized to the weighted number of galaxies brighter than Mj^ s = —21.5 
and Mj = —20.5, divided by the number of groups or clusters as appropriate in the lower two panels. The dashed lines show the (arbitrarily 
renormalized) luminosity function for the whole sample for comparison. In most cases, this curve is hidden by the overlapping luminosity 
functions. 




Fig. 9. — The 67% confidence error ellipses of the J-band Schcchter fits for galaxies divided by environment. On the left are results for 
galaxies in the field (top), group (middle) and cluster (bottom) samples. The heavy ellipses correspond to the Schechter parameters for that 
bin; the lighter ellipses correspond to the Schcchter parameters for the total sample, reproduced in each panel for comparison. On the right 
side of the figure, the samples are divided into EL and NEL galaxies. The best fit values for the emission line galaxies are indicated with 
a square, and the values for the NEL galaxies with a cross. All parameters are for Q = 1, A = 0, and include only k— corrections, with no 
evolution. 




Environment 


Spectral Type 


M* K 


a K s 


M} 


aj 


Total 


all 


-23.48±0.08 


-1.10±0.14 


-22.23±0.07 


-0.96±0.12 




EL 


-23.34±0.16 


-1.46±0.26 


-22.13±0.14 


-1.39±0.19 




NEL 


-23.46±0.12 


-0.81±0.28 


-22.17±0.09 


-0.59±0.15 


Field 


all 


-23.43±0.12 


-1.12±0.21 


-22.18±0.10 


-1.04±0.18 




EL 


-23.28±0.20 


-1.48±0.36 


-22.07±0.17 


-1.54±0.27 




NEL 


-23.38±0.16 


-0.69±0.28 


-22.06±0.10 


-0.44±0.24 


Groups 


all 


-23.58±0.13 


-1.14±0.26 


-22.37±0.17 


-0.88±0.20 




EL 


-23.62±0.44 


-1.64±0.43 


-22.32±0.29 


-1.34±0.34 




NEL 


-23.49±0.26 


-0.81±0.33 


-22.18±0.17 


-0.54±0.27 


Clusters 


all 


-23.81±0.40 


-1.30±0.43 


-22.77±0.34 


-1.30±0.32 




EL 


-23.31±0.70 


-1.18±0.76 


-22.43±0.65 


-1.37±0.56 




NEL 


-23.91±0.52 


-1.28±0.50 


-22.80±0.43 


-1.22±0.39 



Table 1 

We use the cosmological parameters fi = 1, A = 0, h = 1. The tabulated la errors are strongly correlated, 

SHOWN BY THE ERROR ELLIPSES IN FIGURES ^, || AND 
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perfect tracer of the stellar mass; the mass-to- 
light ratio (M/L) can be expected to increase by 
about a factor of two between early- and late-type 
galaxies (e.g., |Bell and de Jong 2001|). To esti- 



mate M/L, constraints on the mean stellar age 
are required. The optical-IR colors of our sample 
are unreliable due to the shallow isophotal limit of 
the LCRS, and the IR colors alone are insensitive 
to stellar population differences. The best indica- 
tor available to us is therefore the -D4000 spectral 
index from the LCRS spectra. We use the solar- 
metallicity models of Fioc & Rocca-Volmerange 
(|1997|) to generate spectra from a given star for- 
mation history. We then compute -D4000 directly 
from the spectra, using the same definition as 
Zabludoff et al. ( |1996j ). The correlation between 
Mj Lj and -D4000 depends only weakly on star for- 
mation history. Figure [11] shows model results 
for a single burst, simple stellar population (SSP) 
and for a galaxy with a constant star formation 
rate, using a Kennicutt ( |1983| ) initial mass func- 
tion. The two models can differ by as much as 
40%, over a narrow range of -D 40 oo, but the rank- 
ing with Mj L is unchanged. -D4000 is sensitive to 
metallicity and reddening, which we neglect, but 
is less sensitive than optical colors. For example, 
an extinction of Ay = 1.5 magnitudes only red- 
dens -D4000 by about 0.2. 

In Figure |12| we show the 67% and 95% error el- 
lipses on the Schechter parameters for the stellar 
mass functions derived using the M/L corrections 
derived from the SSP model. The differences ob- 
served in the luminosity functions are exacerbated 
in the stellar mass functions. Thus, we conclude 
that the stellar mass function is not universal, but 
depends on both galaxy type and the local envi- 
ronment. 

4.3. Implications and Speculations 



Both C hristl ein (|2000D and Zabludoff & 
Mulchaey ( P000|) find that the change in the op- 
tical LF shape with environment is due almost 
exclusively to the NEL galaxies. We confirm that 
the difference in the infrared cluster LF, relative 
to the field, is dominated by this population. It 
is interesting that the faint end slope of the NEL 
galaxies in clusters, as parametrized by a, is simi- 
lar to that of the overall field faint end slope. One 
interpretation for this similarity is that the bulk 
of the cluster population is built up by accreting 
field galaxies with little effect other than the ces- 
sation of star formation, as modeled for example 
by Balogh, Navarro & Morris (|2000|) . However, 
from only the agreement of the LFs, we cannot 



exclude scenarios in which the similar initial pro- 
gentors of faint galaxies in various environments 
began forming stars at different times (the cluster 
ones earlier by virtue of originating near a large 
mass perturbation). In the latter scenario, one 
expects the field LF to evolve to the cluster LF 
in time, as gas is consumed by galaxies and their 
star formation stops. 
In Figure [IB] we show the fraction of NEL galax- 
function of luminosity for the global sam- 
ple. This Figure shows a strong trend that, unfor- 
tunately, is difficult to interpret because metallic- 
ity and extinction also correlate with luminosity 
in a way that makes [Oil] stronger in low lumi- 
nosity galaxies (e.g., Jansen |2001| ). We cannot 
infer the extent to which this relation indicates 
an inherent correlation between galaxy mass and 
instantaneous star formation rate — a correlation 
that would be an important test of galaxy forma- 
tion models. In particular, a variation in the frac- 
tion of currently star forming galaxies with mass 
could indicate that the dominant mode of star for- 
mation (i.e., burst-like or continuous) for a galaxy 
in the field is mass dependent (e.g., Kauffmann, 
Chariot & Balogh, in preparation). 

5. CONCLUSIONS 

We have measured the J and i^-band galaxy 
luminosity functions from the 2MASS, using red- 
shifts from the LCRS, and divided the sample 
based on spectral type and environment. We draw 
the following conclusions: 

• In non-cluster environments, the infrared 
luminosity function of emission line (EL) 
galaxies is steeper than that of non-emission 
line (NEL) galaxies. 

• There is a statistically significant difference 
between the cluster and field galaxy infrared 
luminosity functions, in the sense that M* is 
brighter, and the faint end slope a is steeper, 
in clusters. 

• Using the -D4000 spectral index from the 
LCRS spectra to aid us in converting from 
infrared luminosity to stellar mass, we con- 
clude that the results just described also hold 
for the stellar mass functions and, in fact, the 
differences are more pronounced. 

• The shape of the NEL galaxy stellar mass 
function in clusters is similar to that of the 
global stellar mass function in the field. One 
explanation of this trend is that the cluster is 
mostly built up of field EL galaxies, in which 
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1.2 1.4 1.6 1.8 2 2.2 

D4000 

Fig. 11. — Mass-to-light ratio in the J band, as a function of D400O1 for a Kennicutt (1983) initial mass function from the models of Fioc & 
Rocca-Volmerangc (1997). The solid line is a single burst population, and the dashed line represents a galaxy with constant star formation 
rate. In both models, the galaxy evolves from left to right with time, reaching its maximum -D4000 at 13 Gyr. 




Fig. 12. — 67% and 95% confidence ellipses on the Schechter parameters for various stellar mass functions, using J-band mass-to-light 
ratios determined from the D4000 index, as described in the text. The point with error bars is the result of Cole et al. (2000b) using the 
Kennicutt (1983) initial mass function. The differences in luminosity function shape when divided by spectral type and environment are even 
more pronounced in the stellar mass functions than in the luminosity functions. 
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Fig. 13. — The fraction of non-emission line galaxies as a function of luminosity, as determined by the non-parametric luminosity function 
(points) and the Schechter function fits (solid line). 



star formation has ceased either because (1) 
the galaxies consumed all their gas (either 
through accelerated star formation due to 
the cluster environment, or by forming ear- 
lier due to their proximity to a large mass 
perturbation); or (2) they lost some or all 
their gas through processes like ram pressure 
stripping. 
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